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ABSTRACT 

A new self-sustained model for long-term light curve variations of RX J0513.9— 6951 is proposed based 
on an optically thick wind model of mass-accreting white dwarfs (WDs) . When the mass accretion rate 
to a WD exceeds the critical rate of ~ 1 x 10~^Mq yr~^, optically thick strong winds begin to blow from 
the WD so that a formation of common envelope is avoided. The WD can accrete and burn hydrogen-rich 
matter atop the WD at the critical rate. The excess matter transferred to the WD above the critical 
rate is expelled by winds. This is called the accretion wind evolution. This ejection process, however, 
occurs intermittently because the mass transfer is attenuated by strong winds: the strong winds collide 
with the secondary surface and strip off the very surface layer of the secondary. The matter stripped-off 
is lost from the binary system. Properly formulating this mass stripping effect and the ensuing decay 
of mass transfer rate, we are able to reproduce, in a self-sustained manner, the transition between the 
optical high/X-ray off and optical low/X-ray on states of RX J0513.9-6951. Thus RX J0513. 9-6951 
is the first example of the accretion wind evolution, which is a key evolutionary process in a recently 
developed evolutionary path to Type la supernovae. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: individual 
(RX J0513.9— 6951) — stars: winds, outflows — X-rays: stars 



1. INTRODUCTION 

The Large Magellanic Cloud (LMC) transient supersoft 
X-ray source RX J0513.9-6951 (hereafter RX J0513) was 
discovered by ROSAT (Schaeidt et al. 1993) and has been 
extensively observed (e.g., Reinsch et al. 2000, for recent 
progress). Soon after the X-ray discovery, it was optically 
identified as a binary with an orbital period of 0.76 days 
(Pakull et al. 1993). Its remarkable observational fea- 
tures are summarized as follows: (1) Optical monitoring 
of RX J0513 shows quasi-regular transitions between high 
{V - 16.6) and low {V - 17.4) states (Alcock et al. 1996). 

(2) The duration of the optical high state is ^ 60—150 days 
while the low states last for ^ 40 days (see, e.g.. Fig. 1 of 
Cowley et al. 2002, for recent observational summary). 

(3) Copious supersoft X-rays (~ 30 — 40 eV) were detected 
only in the optical low state (e.g., Schaeidt et al. 1993; 
Reinsch et al. 1996; Southwell et al. 1996). (4) The tran- 
sitions between X-ray off and on are very rapid, in about 
one day or so, while the optical transitions from high to 
low occurs in several days (e.g., Reinsch et al. 2000). The 
transition mechanism between high and low states has not 
been fully elucidated yet, although a few ideas were pro- 
posed so far (see, e.g., Reinsch et al. 2000, for a recent 
progress). 

Very recently, based on an optically thick wind model 
of mass- accreting white dwarfs (WDs), Hachisu & Kato 
(2003) have proposed a new transition mechanism between 
optical high/X-ray off and optical low/X-ray on states. 



They first modeled light curves of the 2000 outburst of 
the recurrent nova CI Aql and then showed that the size 
of a disk around the WD expands largely up to the com- 
panion star or over in a strong wind phase. Including this 
effect of disk expansion in the strong wind phase, they have 
reproduced a sharp ~ 1 mag rise in the optical light curve 
of RX J0513. This is the optical high/X-ray off state, be- 
cause the optically thick winds of WDs certainly obscure 
supersoft X-rays. On the other hand, it drops sharply by a 
magnitude when the wind stops because the size of the disk 
returns to its original size. This is the optical low/X-ray on 
state. However, Hachisu & Kato (2003) gave only a sim- 
ple outline of their mechanism without detailed physical 
descriptions. In this paper, we propose a full description 
of the mechanism and calculate light curves for various 
conditions. We are able to reproduce basic features of the 
observed light curves by this new mechanism. 

It has been frequently discussed that supersoft X-ray 
sources (SSSs) are progenitors of Type la supernovae (SNe 
la) . One of the promising evolutionary models to SNe la 
has been proposed by Hachisu, Kato, & Nomoto (1996). 
They theoretically found that mass accreting WDs blow 
strong winds when the mass accretion to the WD exceeds 
a critical rate, i.e., Afacc > Ma ^ 1 x IO^^A/q yr~^. 
The so-called common envelope evolution does not occur 
(see Hachisu & Kato 2001, for a recent summary). The 
WD accretes hydrogen-rich matter and burns it atop the 
WD at the critical rate of Mcr- The excess matter trans- 
ferred above the critical rate is expelled by winds, i.e.. 
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-^wind = Afacc " -^cr- In such a situation, the WD can 
grow up to near the Chandrasekhar mass (and eventually 
explode as an SN la). This is called the accretion wind 
evolution (Hachisu 2002; Hachisu & Kato 2001) against 
the common envelope evolution (Ibcn & Tutukov 1984; 
Webbink 1984). Thus, the accretion wind evolution is a 
key evolutionary process to produce Type la supernovae. 

Binary systems in the accretion wind evolution are 
predicted to exist but no such systems have been iden- 
tified yet. Here, we have reached the conclusion that 
RX J0513.9— 6951 is the first example just in the accretion 
wind evolution, i.e., in the way to an SN la, the reason of 
which will be discussed below. A brief summary follows: 
The wind ejection process occurs not continuously but in- 
termittently because the mass transfer is attenuated by 
the strong wind itself. 

In §2, an outline of our basic model of RX J0513 is 
briefly introduced and summarized. Our new mechanism 
of transition between the optical high and low states, in 
which the mass transfer rate is attenuated by winds, is 
formulated in §3. A limit cycle model is given in §4 for 
RX J0513 long-term light curve variations. Discussion and 
conclusions follow in §5 and §6, respectively. 

2. BRIEF DESCRIPTION OF A LIMIT CYCLE MODEL 

Long-term variations in the optical light curves of 
RX J0513.9— 6951 are considered as a limit cycle (e.g., 
Reinsch et al. 2000). Before going to detailed descrip- 
tions, we here give a rough sketch on our new mechanism. 
We adopt a binary system consisting of a WD, a disk 
around the WD, and a lobe-filling main-sequence (MS) 
companion, as illustrated in Figure 1. A circular orbit is 
assumed. 

An essential thing is that WDs blow an optically thick, 
strong wind when the mass accretion rate to the WD 
exceeds the critical rate, i.e., Mace > -^cr ^ 1 x 
10~^ Mq yr-^ If the MS companion is 2 - 3 Mq, 
the mass transfer rate from the MS, Mms, is as large as 
a few times 10~^Mq yr~^ or more (e.g., Li & van den 
Heuvel 1997). We distinguish Mms from Mace because 
Mms 7^ Mace in some cases mentioned below. 

Time-evolution of our limit cycle model is schematically 
summarized in Figure 2 (see also Hachisu & Kato 2003). 
A rapid mass accretion to the WD begins at A in Figure 
2. The mass of the WD envelope rapidly increases and the 
WD envelope expands to blow an optically thick, strong 
wind. The strong wind causes a quick rise of V-mag at 
B, because the disk expands to a few to several times the 
previous size, that is, the disk expands from configuration 
(c) to configuration (a) in Figure 1. Because the disk is 
the main source in the optical light due to its irradiation 
effect, this expansion of the disk area causes a magnitude 
rise in the optical light. The massive winds having the 
rate of Mwind ~ 1 x 1O~^M0 yr~^ certainly absorb su- 
persoft X-rays. Thus, this phase is the optical high and 
X-ray off state. As the WD photosphere expands further 
during the rapid mass accretion phase, optical brightness 
increases further from B to C even if the disk shape does 
not change in the wind phase. 

The strong wind hits the MS companion and then strips- 
off its very surface layer. The gas stripped-off is lost from 
the binary. Therefore, the rapid mass accretion is sup- 



pressed by the strong wind itself. The mass-stripping rate 
is roughly proportional to the wind mass loss rate. As 
the wind mass loss rate increases, the mass-stripping rate 
eventually overcomes the mass transfer rate. This happens 
at C. However, the accretion disk still continues to supply 
mass to the WD at least during the viscous timescale of 
the accretion disk, fyis- Then, mass accretion to the WD 
stops at C about tyis after the mass transfer from the MS 
stops. 

The mass of the WD envelope is now gradually ckx'reas- 
ing due mainly to wind mass loss. The photosphere of the 
WD is also gradually shrinking. This makes a gradual de- 
cay of the optical light curve even if the disk shape is the 
same. The wind eventually stops at D. We expect copious 
supersoft X-rays after the wind stops. The disk shrinks 
to a normal size at E, that is, from configuration (a) to 
configuration (b) in Figure 1. This causes a sharp ~ 1 
mag drop in the optical light. A rapid mass transfer from 
the MS to the disk resumes at E, and then the disk edge 
flares up at F, that is, configuration (b) to configuration 
(c) in Figure 1. It takes several dynamical timescales to 
change from E to F. This phase corresponds to the optical 
low and X-ray on state. 

Because mass accretes to the WD via an accretion disk, 
this rapid mass accretion flow reaches the WD surface at 
G(=A) about tyis after the rapid mass transfer from the 
MS resumes. Thus, the cycle is repeated. 

3. THE BINARY MODEL 

3.1. Basic Parameters of the Binary 

A relatively massive white dwarf is suggested because 
of a rapid contraction timescale of the WD photosphere 
during the transition from optical high to low states (e.g., 
Southwell et al. 1996). We have adopted Mwd = 1.3 
after Hachisu & Kato (2003) but later examined other 
cases of the WD mass, i.e., Mwd = 1-35 Mq, 1.2 Mq, 
and 1.1 Mq. Since the LMC metallicity is reported to 
be about a third of the solar metallicity, we assume the 
metallicity of Z = 0.004 and hydrogen content of X = 0.7 
for the envelope of the WD. The mass ratio is assumed to 
be g = Mms /Mwd = 2 because mass transfer occurs on a 
thermal timescale and a rate of Mms ^ 1O~^M0 yr~^ is 
suggested by Southwell et al. (1996). Then, the compan- 
ion mass is Mms = 2.6Mq for Mwd = 1-3 Mq. The non- 
irradiated surface temperature of the lobe-filling 2.6Mq 
MS companion is 13,000 K for the metallicity of Z = 0.004 
(e.g., Fagotto et al. 1994). 

For the orbital period of Porb = 0.763 days (see below), 
the separation is a = 5.53 Rq, the radii of the effective 
Roche lobes arc Rl = 1.77 Rq for the primary (WD) and 
R2 = R2 = 2.43 Rq for the secondary (MS) component. 
The radius of a 2.6 Mq zero-age main-sequence (ZAMS) 
star is -Rzams = 1-5 Rq for Z = 0.004, so that the sec- 
ondary has already evolved to expand but still remains 
at the main-sequence (before exhaustion of central hydro- 
gen). 

The ephemeris of RX J0513 has been revised by Cowley 
et al. (2002) as 

i(JD) = 2, 448, 858.099 + 0.7629434 x E (1) 

at minima of the optical orbital modulations. We use this 
revised ephemeris with the MS companion being in front 
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of the WD at minima (sec also Fig. 5 of Hutchings et al. 
2002). The iiichnation angle is not known but suggested 
to be low. Here we adopt i = 10°. Calculated light curves 
are not so different among i = 0° — 35° as shown later. 

3.2. Time Evolution of the WD Envelope 

Time evolution of the photospheric radius and temper- 
ature of the WD are calculated from Kato & Hachisu's 
(1994) optically thick wind solutions after Hachisu & Kato 
(2001), i.e., 



at 



Mnuc - -^wind, 



(2) 



where AMcnv is the hydrogen-rich envelope mass of the 
WD, Mace the mass accretion rate to the WD, M„uc the 
decreasing rate by nuclear burning, and Mwind the wind 
mass loss rate of the WD. The nuclear burning rate Mnuc, 
wind mass loss rate Mwind, photospheric radius i?ph, pho- 
tospheric temperature Tph, and luminosity iph as well as 
the photospheric wind velocity fph, of the WD are all given 
as a function of the envelope mass (AMe„v) in the wind 
solutions (Kato & Hachisu 1994; Hachisu & Kato 2001), 
i.e., 

Mnue = /l(AMenv, X, Z, Mwd), (3) 



Mwind = /2(AMenv, X, Z, Mwd), 

i?ph = /3(AMe„v, X, Z, Mwd), 
2ph = /4(AAfonv, X, Z, Mwd), 
iph = fsi^M^nv, X, Z, Mwd), 
Vph = /6(AMe„v, X, Z, Mwd). 



(4) 
(5) 
(6) 
(7) 
(8) 



Given the mass accretion rate Afacc and the initial enve- 
lope mass, we are able to follow the time evolutions of the 
envelope mass and then the radius and temperature (and 
also the luminosity) of the WD by integrating equation 
(2). 

3.3. Optically Thick Strong Winds 

Mass transfer in RX J0513 probably proceeds in a 
Kelvin- Helmholtz timescale of the MS companion (e.g., 
van den Heuvel et al. 1992), so that it is as large as 
Mms 1 X 10~^Mq yr"'^ or more. When the mass trans- 
fer rate, Mms, exceeds the critical rate, i.e., 

Mms > M„ « 0.75 x lO"*^ (^^^ - 0.4^ Mg yr'^, (9) 

the mass-accreting WD expands to blow optically thick 
winds (e.g., Hachisu et al. 1996) and the mass loss rate 
approximately given by 

Mwind ~ Mms - M„, (10) 

when all the transferred matter accretes to the WD, i.e.. 
Mace = Mms- Equation (9) is the same form as that for 
the solar metallicity {Z = 0.02), but valid even for a lower 
metallicity of Z = 0.004 (e.g., Hachisu & Kato 2001). 



3.4. Mechanism of Intermittent Mass Accretion 

When strong winds collide with the surface of the com- 
panion, the surface is shock-heated and ablated in the 
wind. We estimate the shock-heating by assuming that 
the velocity component normal to the surface is dissipated 
by shock and its kinetic energy is converted into the ther- 
mal energy of the surface layer. The very surface layer of 
the envelope expands to be easily ablated in the wind. We 
regard that gas is ablated and lost from L3 point (L3 is 
the outer Lagrange point near the MS companion) when 
the gas gets the same amount of thermal energy as the 
difference of the Roche potentials between the MS surface 
and L3 point. Then the mass stripping rate is given by 

(?!'L3 - </'Ms) • Mstrip = ^■y^Mwind " VeS ' 9{q), (H) 

where M = Mwd + Mms, a is the separation of the binary, 
(t>MS the Roche potential at the MS surface, (/)l3 means the 
Roche potential at L3 point near the MS companion, both 
of which are normalized by GM/a, ries the efficiency of 
conversion from kinetic energy to thermal energy by shock, 
g{q) is the geometrical factor of the MS surface hit by the 
wind including the dissipation effect (normal component 
of the wind velocity), and g{q) is only a function of the 
mass ratio q = Mms /Mwd (see Hachisu, Kato, & Nomoto 
1999a, for more details). Here we modified equation (21) 
of Hachisu et al. (1999a) to include the effect of Roche 
lobe overflow from L3 point. Then the stripping rate is 
estimated as 



M 



strip 



ciM, 



wind. 



where 



^ Vcs ■ gjq) f v'^a \ 

~ 0L3-</'MS \2GMJ 



(12) 



(13) 



The efficiency of conversion is assumed to be rjes = 1. 
Then, we have 



ci w 0.1 



400 km s" 



10, 



(14) 



for Mms < Mgtrip 
for Mms > Mstrip 



(15) 



for the wind velocity of Vy,ind = 3800 km s ^ (e.g.. Cramp- 
ton et al. 1996; Hutchings et al. 2002; Southwell et 

al. 1996), a = 5.5i?0, M = 3.9M0, 0l3 - 0ms = 0.3, 
and g{q) — 0.025 for the mass ratio q = Mms/Mwd = 2 
(Hachisu et al. 1999a). 

In such a wind phase, the net mass accretion rate to the 
WD, Mace, is modified as 

Mace = |'^-("^-^)'- 

I Mms - Mstrip, 
as illustrated in Figure 3. Here we adopt a small value 
of e = 1 X 10"'' Mq yr^^ when Mms < Mgtrip because 
the mass accretion to the WD does not stop abruptly but 
probably continues at least for a draining time of the ac- 
cretion disk. We do not know the exact draining time of 
the accretion disk after the mass transfer from the MS 
stops. Alternatively, we just assume a small rate of the 
mass accretion e to mimic the draining of the disk during 

that Mms < Mtrip- 

To know when the rapid mass accretion resumes, we 
monitor the upper level of flow in Figure 3 as 

^Mflow = Mms - Mstrip - Mace, (16) 
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with the initial value of Mflow = 0. During that Mgtrip > 
MmS) -^flow decreases to a large negative value as illus- 
trated in Figure 3b. Then the mass transfer from the MS 
stops or its rate drops to e. Once the mass accretion rate 
drops below the critical rate, i.e., Mace < Mcr, the wind 
gradually weakens and eventually stops. The stripping ef- 
fect vanishes (Mgtrip = 0) and then the level of flow begins 
to rise. We start again the rapid mass accretion when the 
level reaches Mgow = as illustrated in Figure 3a. 

It should be noticed that equation (2) is on the WD but 
equations (11) — (16) are on the MS. Therefore, the mass 
accretion rate, Mace, may not be the same on each side at 
the specific time, t, because it takes a viscous timescale of 
the accretion disk for gas to reach the WD surface from 
the MS surface. Here we adopt a viscous timescale of 

tvis = -^ = 5-~40 — — days, 17) 

V avis CsH \ 0.1 / 

where i/ is the viscosity, ayis the a parameter of Shakura 
& Sunyaev (1973), Cs the sound speed, and H the vertical 
height of the accretion disk. We adopt i?disk = 1-4 Rq 
and ff/Edisk ~ 0.1 for Tdisk ~ 30000 - 50000 K and 
Cs ^ 30km s~^. Then, we have 



Macc(t) =Macc(i-iv 



WD 



MS 



(18) 



3.5. Irradiation Effects of Disk 

The irradiation effects both of the disk and of the MS 
companion play an essential role in the light curve mod- 
eling. For the basic structure of the accretion disk, we 
assume an axi-symmetric structure with the size and thick- 
ness of 

i?disk = aRl (19) 

and 

(20) 



h = (3R. 



disk 



isk 



where i?disk is the outer edge of the accretion disk, i?J is 
the effective radius of the inner critical Roche lobe for the 
WD component, h is the height of the surface from the 
equatorial plane, and w is the distance on the equatorial 
plane from the center of the WD. Such binary models are 
illustrated in Figure 1. We assume that the photosphere 
of the WD, the MS companion, and the accretion disk 
emit photons as a blackbody at a local temperature which 
varies with position. 

We adopt v = 1 during the strong wind phases, i.e., 
configuration (a) in Figure 1 and just after the wind stops 
(but before a rapid mass transfer from the MS resumes), 
i.e., configuration (b) in Figure 1, but v = 2\n the no wind 
case, i.e., configuration (c) in Figure 1. This n7-square law 
mimics the effect of elevated rim of the accretion disk by 
spray as discussed by Schandl, Meyer-Hofmeister, & Meyer 
(1997) for luminous supersoft X-ray sources. On the other 
hand, we simply assume v = \ va the wind phase, be- 
cause such an elevated rim is blown in the wind during 
a strong wind phase. The efficiency of irradiation is as- 
sumed to be 0.5 both for the disk and for the compan- 
ion (see also Schandl, Meyer-Hofmeister, & Meyer 1997; 
Meyer-Hofmeistcr, Schandl, & Meyer 1997). 

The size and thickness of disk, a and /3, are also the im- 
portant parameters that determine the brightness of disk 



by irradiation effects. Hachisu & Kato (2003) have de- 
termined a and /3 from the light curve fitting of CI Aql 
in outburst. The disk expands up to the companion star 
as shown in Figure la, i.e., a ~ 2 — 3, because the disk 
surface is blown in the wind. It shrinks to a normal size 
of a = 0.7 — 0.8 after the wind stops. These two param- 
eters of (a,/3) are determined as (a, /3) = (3.0,0.05) — > 
(0.8,0.05) (0.8,0.3), during the wind phase, just after 
the wind stops, and then the rapid mass transfer resumes, 
respectively (see Hachisu & Kato 2003, for more details). 
The shrink from a large, extended size to a normal size is 
important in the modeling of a sharp ~ 1 mag drop in the 
optical light curve. Thus, a large expansion of the disk in 
the strong wind phase has been first confirmed in the 2000 
outburst of the recurrent nova CI Aql. 

We find essentially the same sets of (a, (3) as in Hachisu 
& Kato (2003) is valid for RX J0513. We may understand 
such a large expansion of the disk beyond the Roche lobe 
as follows: Observationally very fast winds have been re- 
ported, for example, 3800 km s^^ in RX J0513 (see, e.g., 
Crampton et al. 1996; Hutchings et al. 2002; Southwell 
et al. 1996). The wind mass loss rate is calculated to be 
10~® — 1O~'*M0 yr"-'^ in such a strong wind phase (e.g., 
Hachisu & Kato 2001). There is a large velocity differ- 
ence between the wind and the disk surface. It certainly 
drives the Kelvin-Helmholtz instability at the interface, 
so that the very surface layer of the disk is dragged away 
outward with the velocity at least several hundred km s^"'^, 
i.e., almost a free stream going outward. This surface free 
stream is initially optically thick at the original disk region 
but becomes optically thin somewhere outside because of 
geometrical dilution effect. Only the optically thick region 
of disk stream can contribute to the optical light by irra- 
diation, so that we regard the radius of transition region 
between optically thick and thin as the edge of the disk. 
In this sense, the disk edge as shown in Figure la is not 
the boundary of matter but just a transition from optically 
thick to thin of the disk surface flow. 

It should be noted that a high density part of the disk 
is hardly changed by this relatively tenuous surface free 
stream and resides in within its Roche lobe because the 
internal density of the disk is much denser than that of 
the WD wind. In the case of Figure 4, for example, the 
wind mass loss rate is about 1 x 1{)~^Mq yr~^ and its ve- 
locity is ~ 4000 km s~^, so that the density of the wind 
is estimated from the continuity (Af„ind = 47rr^yOw) to be 
about 10~^^g cm^ at the distance of 1 Rq from the center 
of the WD. On the other hand, the density of the stan- 
dard accretion disk is about 10~^g cm'^ at the same ra- 
dius. Here, we assume the standard accretion disk model 
(Shakura & Sunyaev 1973) with the mass accretion rate 
of ~ 5 X IO-'^Mq yr-i. 

It should be also addressed that optically thick winds are 
accelerated deep inside the photosphere (r 3> 1) but the 
wind itself becomes optically thin (r <C 1) far above the 
WD photosphere as the definition of photosphere. There- 
fore, we do not expect any large optical contribution from 
the wind itself far outside the WD photosphere (see, e.g., 
Meyer-Hofmeistcir ct al. 1997, for an estimation of small 
wind contributions to the optical light). 

The numerical methods for obtaining light curves, in- 
cluding the above irradiation effects, have been described 
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in more detail in Hachisu & Kato (2001). 

4. LIMIT CYCLE MODEL OF LONG-TERM VARIATIONS 

4.1. A template light curve 

Assuming Mms = 5 x 1O~''M0 yr~\ ci = 10.0, and 
tvis = 20.5 days, we have calculated time evolution of our 
binary system with the initial envelope mass of AMcnv = 
5 X 10~^AfQ. The obtained V-light curve and color (second 
panel), accretion rate (Afacc on the WD) and wind mass 
loss rate in units of lO^^M© yr"-*^ (third panel), radius 
(-R0) and temperature (eV) of the WD photosphere (bot- 
tom panel) together with the MACHO observation (top 
panel) are shown in Figure 4. Direct fitting to the obser- 
vational light curves gives us an apparent distance mod- 
ulus of (m — M)v = 18.7. If we use the color excess of 
E(B -V) = 0.13 (Gansicke et al. 1998), the distance to 
RX J0513 is (m — M)o = 18.3, being consistent with the 
relatively short distance to the LMC. In this paper, we use 
(to - M)v = 18.7. 

At the beginning of time-evolution calculations, the 
mass of the WD envelope is small so that the WD does 
not blow a wind yet. Its envelope mass increases because 
of Mms = Mace > Mnuc (Mwind = 0). When the envelope 
mass reaches AMenv = 6.5 x 1O~^M0, the WD photo- 
sphere expands to i?ph = 0.087 Rq and begins to blow a 
wind. The disk size expands from (a, /3) = (0.8,0.2) to 
(3.0, 0.05) in dynamical timescales. Here we adopt a few 
days as the transition time, because the strong winds eas- 
ily blow the very surface layer of the disk in the wind. The 
wind mass loss rate is about Mwind = 1 x 1O~^M0 yr~^ 
only one day after the wind starts, which is large enough 
to completely obscure soft X-rays (see, e.g., discussion of 
Southwell et al. 1996). The occurrence of optically thick 
winds can naturally explain the very rapid emergence and 
decay (< 1 — 2 days) of supersoft X-rays accompanied 
by the transition between the optical high and low states 
(Reinsch et al. 2000). We discuss again in more detail 
this transition timescale of soft X-rays bellow in §5. 

The envelope mass continues to increase and the wind 
mass loss rate also grows. The WD photosphere also 
continues to expand and this causes a gradual increase 
in the optical light curve. The mass transfer from the 
MS is suppressed when the wind mass loss rate reaches 
Mwind = 0.5 X 1O"'^M0 yr-\ because M^trip = lOMwind = 
Mms — 5.0 x 10~^ Mq yr~^. However, the mass accretion 
to the WD still continues at a high rate and the wind mass 
loss rate reaches M^ind = 1-0 x 1O~^M0 yr~^ as shown in 
Figure 4, because the mass in the accretion disk is drained 
to the WD in another 20.5 days of t^is time delay. The 
mass accretion to the WD eventually stops and the mass 
of the WD envelope begins to decrease due to wind mass 
loss and nuclear burning. The WD photosphere gradually 
shrinks from i?wD,ph — 0.24 Rq to 0.087 Rq during the 
wind phase. This causes the gradual decrease in the opti- 
cal light curve until the end of high state. When the mass 
of the WD envelope decreases to AMenv = 6.5 x 1O~^M0, 
the wind stops. 

The WD photosphere begins to shrink quickly after the 
wind stops. It takes 34 days to collapse from 0.087 Rq to 
0.03 Rq. The photospheric temperature increases quickly 
from 30 eV to 40 eV during this period, which is very con- 



sistent with the blackbody temperature of supersoft X-rays 
(Schaeidt et al. 1993). When the wind stops, the wind 
mass loss rate decreases from M^ind = 1-0 x 1O~^M0 yr"-"^ 
to zero in a day. This is again very consistent with a 
sharp emergence of supersoft X-rays (about a day or so, 
see Reinsch et al. 2000). 

In Figure 3, during the strong wind phase the schematic 
mass flow goes to much below the level of full. After the 
wind stops, it gradually goes up and eventually reaches 
the level of full. It takes 13.5 days in Figure 4. Then, 
the mass accretion to the WD resumes 20.5 days after it 
reached the level of full, i.e., tvis delay time. It takes more 
6 days that the WD envelope expands to blow a wind 
again. The total duration of the low state amounts to 
13.5 -|- 20.5 -I- 6 = 40 days. Thus, the system repeats the 
cycle mentioned above. 

The modeled system reaches a limit cycle after one cy- 
cle of high and low states. The long-term light curve 
modulations are well reproduced. The time-averaged mass 
transfer rate is 1.1 x 1O~^M0 yr~^ and 35% of the trans- 
ferred matter is blown in the wind, so that the MS loses 
its mass at the average rate of 3.9 x 10~^Mq yr~^ by 
stripping effects. These values satisfy the conditions of 

Mms = Mace + Mgtrip and M^trip = 10 M^ind in averaged 
values. 

It should be noted that our ivis adopted here is 20.5 
days, which is tuned to reproduce the duration of optical 
low states, i.e., ^ls = 40 days, because the exact number 
of tvis is still theoretically not known to us. The duration 
of low states shortens to iLS — 25 days when tvis = 18 
days. Optical low states disappear for tvis ^ 14 days. On 
the other hand, it becomes longer ths = 50 days when 
tvis = 23 days. We summarize these results in Table 1. 
The effect of tvis on the long-term light curves will be dis- 
cussed more below in relation to other system parameters. 

Thus, we are able to reproduce the basic observational 
features (1) — (4) of RX J0513 summarized in §1. Fur- 
thermore, our model naturally explains two other observa- 
tional features: (5) the WD photosphere continues to ex- 
pand during the rapid mass accretion phase and still goes 
on until ~ 20 days after the wind started. This makes 
a bump of ~ 0.1 mag at the head of optical high states. 
(6) The WD photosphere gradually shrinks after the rapid 
mass accretion stops. This also make a gradual decay of 
~ 0.2 mag until the end of the wind phase (until the end 
of the optical high state). 

4.2. Mass transfer rate of MS 

The mass transfer rate of the MS companion depends 
on its evolutionary stage and binary parameters (e.g., Li 
& van den Heuvel 1997; Langer et al. 2000). Since we 
do not know the exact number of the mass transfer rate in 
RX J0513, we have examined various values of the original 
mass transfer rate Mms, i-e., 0.6, 0.7, 0.8, 0.9, 1.0, 2.0, 3.0, 
5.0 (the template), 7.0, 9.0, 10.0, 15.0, and 20.0 in units 
of 1O~^M0 yr~^, as summarized in Table 1. The duration 
of optical high states becomes longer for the higher mass 
transfer rates of Mms- On the other hand, high/low states 
almost become a 50% duty cycle as the mass transfer rate 
approaches the critical value, i.e., Mms Mcr. Some of 
the calculated light curves are shown in Figure 5. 

A remarkable result is that the durations of low states 
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are kept to be ~ 40 days for a wide range of the mass 
transfer rates, i.e., 7 x 1O~'^M0 yr~^ < Mus $ 6 x 
10~^Mq yr"-'^. The duration of high states varies largely 
from 40 days to 120 days for the same range of the mass 
transfer rates. This is just the result that we want to ex- 
plain theoretically, because it is an observational feature 
as shown in Figure 4 (top panel). The low state gradually 
shortens and eventually disappears as the mass transfer 
rate increases to Mms ^ 1-5 x 1O~^M0 yr~^ for a set of 
tyis = 20.5 days and Ci = 10.0. The luminosity of high 
states is brighter for higher mass transfer rates of Mms- 

4.3. White dwarf mass 

We have also examined various WD masses, i.e., Mwd = 
1.35, 1.2, and 1.1 Mq, with the mass ratio being kept 
q = Mms /Mwd = 2. The results are shown in Figures 
6— 8. In general, more massive white dwarfs have a larger 
wind mass loss rate, so that the duration of the wind phase 
becomes much shorter. Because we do not know the ex- 
act number of tvis, our guided principle is that we tune 
tvis to satisfy the duration of the low state is 40 days, i.e., 
t-Ls = 40 days. For the case of Mwd = 1-35 Mq shown 
in Figure 6, we have tuned the viscous timescale Uis at 
the mass transfer rate of Mms = 5 x IO^'^Mq yr^^. The 
tuned viscous timescale is tvis = 15.0 days. The durations 
of high/low states are similar to those for Mwd = 1-3 Mq 
in Figure 5 except for the relatively short durations. The 
WD mass in RX J0513 may be less massive than 1.35 Mq, 
because the gradual decay of the light curve in high states 
seems to be too fast to be compatible with the observation. 

For the case of Mwd = 1-2 Mq shown in Figure 7, 
we have tuned the viscous timescale at the mass trans- 
fer rate of Mms = 3 x IO^^Mq yi-\ i.e., Uis = 32.0 
days. The durations of high states are longer than those for 
Mwd = 1-3 Mq in Figure 5. The WD mass in RX J0513 
may be somewhere between 1.2 Mq and 1.3 Mq, because 
the decay of the light curve in high states seems to be 
slightly slower than, but in reasonable agreement with, 
the observation. 

For the case of Mwd = 1.1 Mq shown in Figure 8, we 
have tuned the viscous timescale at the mass transfer rate 
of Mms = 2 x 10~^ Mq yr-\ i.e., = 45.0 days. The 
durations of high states are much longer than those for 
Mwd = 1.3 Mq in Figure 5. The WD mass of 1.1 Mq 
may not be the case, because the decay of the light curve 
in high states is too slow to be compatible with the obser- 
vation. 

4.4. Mass stripping effect 

It is a difficult work to accurately estimate the effect 
of mass-stripping, that is, to determine the coefficient of 
Ci in equation (12). Therefore, we calculate two other 
cases of ci, i.e., ci = 5 and ci = 1.5 and examine how 
the light curves depend on the coefficient. It is obvious 
that strong modulation of the mass transfer rate never oc- 
curs for the case of ci < 1 because there is no case of 
Mms < Mstrip = ciMwi„d(< Mwi„d < Mace < Mms)- Our 
calculated models and their parameters are listed in Table 
1. 

For the case of ci =5 shown in Figure 9, the over- 
all features of light curves are essentially similar to the 
case of ci = 10 shown in Figure 5. Main differences are 



that (1) the tuned viscous timescale is <vis = 36.0 days at 
Mms = 2 X 10~^Mq yr~^, much longer than tyis = 20.5 
days for ci = 10, and (2) bumps at the head of high states 
continue for a much longer period and may not be com- 
patible with the observation. 

For the case of ci = 1.5 shown in Figure 10, a new 
thing is the appearance of no modulation in mass trans- 
fer for a range of 0.67 x IO'^'Mq yr^^ < Mms < 1-5 x 
10~^Mq yr~^. The tuned viscous timescale is <vis = 62.0 
days at Mms = 2 x 10~^Mq yr~^. Bumps at the head of 
high states continue for about 40 days. Both the viscous 
timescale and the duration of bumps are much longer than 
those for Ci = 10 and ci = 5. The long duration of bumps 
is not compatible with the observation. 

An important result in this subsection is that modu- 
lation in mass transfer occurs and optical high/low states 
appear even for relatively low efficiencies of mass-stripping 
on the MS, i.e., for ci — 1.5, if the viscous timescale is rel- 
atively long, i.e., tyis = 62 days. Such a low case of ci may 
correspond to a low velocity case of winds as understood 
from equation (14) and a relatively long viscous timescale 
is realized for binaries with a longer orbital period, for ex- 
ample, longer than a day, as derived from equation (17). 

5. DISCUSSION 

It has been argued that only a viable mechanism for 

the transition between the optical high/X-ray off and op- 
tical low/X-ray on states is the contraction/expansion of 
the WD photosphere resulting from the variation of the 
mass accretion rate to the WD (e.g., Meyer-Hofmeister et 
al. 1997; Pakull et al. 1993; Reinsch et al. 1996, 2000; 
Southwell et al. 1996), although a few other mechanisms 
have ever been proposed (see, e.g., discussion of Southwell 
et al. 1996). Therefore, we first examine whether or not 
this mechanism can explain all the observational features. 

Reinsch et al. (2000) detected a very sharp rise (< 3 
days) and a steep decline (< 2 days) in supersoft X-rays by 
a factor of > 100. This requires a contraction/expansion 
of the WD photosphere by a factor of > 7 and an in- 
crease/decrease of the photospheric temperature by a fac- 
tor of > 3 (Reinsch et al. 2000). We have estimated the 
timescale of contracting/expanding WD envelopes based 
on our wind and static solutions. The WD photosphere 
gradually shrinks during the wind phase while it begins to 
contract most rapidly just after the wind stops (see the 
bottom panel of Fig. 4). Even for this most rapid shrink 
phase, it takes about three Kelvin-Helmholtz timescale to 
shrink by a factor of four, i.e., 

icontract ~ 3tkH, (21) 

from our model calculations. This Kelvin-Helmholtz 
timescale is defined by 

^-^env, thermal /rio\ 
TKH = ^ , (22) 

where AE^nv, thermal is the thermal energy of the WD en- 
velope, and L the luminosity of the WD, just when the 
wind stops. These contraction timescales just after the 
wind stops are plotted against various WD masses in Fig- 
ure 11. It is obvious that the rapid (< 2 days) change in 
the soft X-ray flux cannot be explained by the contraction 
of the WD photosphere even if we assume the extreme case 
of Mwd = 1-377 Mq (at least, for Z = 0.004). 
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For the expansion case of the WD photosphere dur- 
ing the rapid mass accretion phase, our steady-state wind 
model does not correctly include the compression term by 
accretion so that the estimated timescale is not so ac- 
curate. In this sense, our numerical results cannot give 
a strong constraint on the timescale of photospheric ex- 
pansion. It should be pointed out, however, that usual 
time-dependent Henyey methods fail to be converged in a 
rapidly expanding photosphere of WDs and are now un- 
able to involve wind mass loss. 

Reinsch et al. (2000) suggested a characteristic de- 
cay timescale of 34 days from the exponential decay of 
supersoft X-rays in an optical low state. In our model, 
there are two possibilities. One is the timescale of expan- 
sion/contraction of the WD photosphere and the other is 
the viscous timescale of the accretion disk, tvis- If this 
34 days timescale is related to the expansion/contraction 
timescale of a WD envelope, its WD mass is somewhere be- 
tween 1.2 Mq and 1.3 M© for the mctallicity oi Z = 0.004 
as shown in Figure 11. This mass range of the WD is 
consistent with our long-term light curve modulations. 
On the other hand, if it indicates that t^is = 34 days, 
the ranges of the parameters are also suggested to be 
MwD = 1-2 - 1.3 Mq and ci = 5 - 10 from Table 1. 
These are also consistent with our long-term light curve 
modulations. 

Southwell et al. (1996) suggested that star spots on 
the secondary surface cover the LI point, resulting in a de- 
creased mass transfer rate from the secondary. This mech- 
anism has been originally suggested by Livio & Pringle 
(1994) for the VY Scl stars. In RX J0513, however, it is 
very unlikely that the MS companion develops star spots 
on its surface, because its mass is suggested to be 2— 3 Mq, 
not having a deep surface convection. 

Very recently, Hutchings et al. (2002) reported radial 
velocity modulations of 117 ± 40 and 54 ± 10 km 
for the broad O VI emission and Lyman absorption, re- 
spectively. If these trace the motions of the hot and 
cool components, respectively, the estimated mass ratio 
is q = Mms/Mwd = 2.2 ± 1.2, being consistent with our 
assumption of g = 2. The corresponding inclination an- 
gle is z = 28° ± 10° for our case of Mwd = 1.3 Mq and 
.^MS = 2.6 Mq, which is a bit larger than our assumption 
of i = 10°. Then, we have examined how the long-term 
light curve depends on the inclination. Additional four 
cases of the inclination angle, i.e., z = 0°, 20°, 28°, and 
35°, are calculated. The results are also shown in Fig- 
ure 12, which are not so much different in shape from the 
case of z = 10°, but the brightness in optical high states 
decreases by a factor of cosz. 

Cowley et al. (2002) have reanalyzed the MACHO data 
for eight years and determined a new orbital period and a 
new ephemeris. They found that the orbital modulations 
of the optical light curve are different in shape between the 
optical high and low states. The orbital modulation dip 
is much broader in the high state than in the low state. 
Assuming an inclination angle of i = 28°, we have cal- 
culated two orbital modulations each for the optical high 
state and for the optical low state as shown in Figure 13. 
The observational broad and shallow orbital modulations 
in the optical high state are reasonably reproduced as well 
as their full extents (amplitudes) of AV 0.06 mag. On 



the other hand, the full extents of the orbital modulations 
become larger (deeper) to AV ~ 0.09 mag in the optical 
low state, being consistent with the observation. 

6. CONCLUSIONS 

A new self-sustained model for the long-term light curve 
variations of RX J0513.9— 6951 is formulated based on our 
optically thick wind model of mass-accreting white dwarfs. 
We also show that RX J0513.9— 6951 is the first beautiful 
example of the accretion wind evolution, which is a key 
evolutionary process to Type la supernovae (e.g., Hachisu 
2002; Hachisu & Kato 2001) 

We have already shown that, when the mass accretion 
rate to a WD exceeds the critical rate M„ w 0.75 x 
1O-6(Mwd/M0 - 0.4)Mq yr^i, optically thick strong 
winds begin to blow from the WD. As a result, a forma- 
tion of common envelope is avoided. In such a situation, 
the WD accretes and burns hydrogen atop the WD at the 
critical rate. The excess matter transferred to the WD 
above the critical rate is expelled by winds (Hachisu et al. 
1996, 1999a,b). This is called the accretion wind evolution 
(Hachisu & Kato 2001). Unfortunately, no such systems 
have ever been discovered. Here, we report that the LMC 
supersoft X-ray source RX J0513.9— 6951 is the first ex- 
ample of such an object, i.e., just in the accretion wind 
evolution. 

We summarize our main results as follows: 

1. Mass ejection in the accretion wind evolution does oc- 
cur not continuously but intermittently because the mass 
transfer is attenuated by strong winds. The strong winds 
collide with the secondary and strip off the very surface 
layer of the secondary. Properly formulating this mass 
stripping effect and the ensuing decay of mass transfer 
rate, we obtain, in a self-sustained manner, on/off of mass 
transfer from the companion star. 

2. A light curve modeling on the 2000 outburst of the 

recurrent nova CI Aql showed that the size of the disk 
around the white dwarf extends widely up to the compan- 
ion star or over only during the strong wind phase (Hachisu 
& Kato 2003). Including this large extension of a disk in a 
strong wind phase, we are able to reproduce the transition 
between the optical high and low states of the LMC super- 
soft X-ray source RX J0513.9— 6951: When a strong wind 
occurs, the very surface layer of the disk is blown in the 
wind and the disk expands up to the companion or over. 
This large extension of the irradiation area makes a sharp 
~ 1 mag rise in the optical light curve. On the other hand, 
it drops sharply by a magnitude when the wind stops be- 
cause the size of the disk returns to its original size. 

3. In optically thick winds, the mass loss rate rises 
to Mwind ?J 1 X 1O~*M0 yr^^ one day after the wind 
starts. This can obscure completely soft X-rays. On the 
other hand, the wind mass loss rate drops from Mwind ~ 
1 X 1O~^M0 yr~^ to zero in within a day, which en- 
ables a short (< 1 day) timescale of soft X-ray emergence. 
Thus, we are also able to reproduce a very rapid emer- 
gence/decay of supersoft X-rays in within a day or so. 

4. We abstract two important system parameters. One 
is the viscous timescale of the accretion disk, t^is, and 
the other is the efficiency of mass-stripping by winds, ci, 
defined by Mgtrip = ci Mwind- Both of them should be 
determined self-consistently but are not accurately deter- 
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mined at least at the present time like the a-parameter 
of Shakura & Sunyaev (1973). We are able to repro- 
duce the optical high/low states for reasonable ranges of 
the system parameters, i.e., 20 days < t^is < 60 days and 
1-5 < ci < 10. 

5. The original mass transfer rate from the main- 
sequence companion is probably constant in time. Even if 
this original mass transfer rate varies from several times 
10~^ Mq yr~^ to several times 1O~^M0 yr~^, by a factor of 
ten, the duration of optical low states hardly changes, i.e., 
stays around ~ 40 days, for a fixed set of (iyisi ci). On the 
other hand, the duration of optical high states varies from 
~ 40 days to 120 days for the same range of the original 
mass transfer rate. These durations of optical high/low 
states arc the prominent features of the MACHO observa- 
tions (Alcock et al. 1996; Cowley et al. 2002). 

6. The white dwarf photosphere expands and blows a wind 
during the rapid mass accretion phase. The photosphcric 
expansion still goes on for a while after the wind suppresses 
the mass transfer, because the accretion disk supplies mass 
to the WD for a draining time (viscous timescale). This 
makes a bump of ~ 0.1 mag at the head of the optical high 
state. The duration of this bump may determine the value 
of tyis- The WD photosphere gradually shrinks after the 
rapid mass accretion stops. This makes a gradual decay 
of ~ 0.2 mag in the optical light curve until the end of the 
wind phase (until the end of the optical high state). These 
two light curve features clearly appear in the MACHO ob- 
servation (Alcock et al. 1996). 



7. Comparison of the modeled light curves with the 
observation implies that the mass of the white dwarf 
is about 1.2 — 1.3Mq, the original mass transfer rate 
(2 - 3) X 1O~^M0 yr~^, the viscous timescale 20 - 40 days, 
and the efficiency of stripping by winds ci = 5 — 10. 

8. The on/off timescales of supersoft X-rays are as short 
as one or two days (Reinsch et al. 2000). Only the con- 
traction/expansion of the white dwarf photosphere cannot 
explain the rapid emergence/decay of supersoft X-rays be- 
cause the contraction/expansion timescales of the white 
dwarf photosphere are too long to reproduce such a rapid 
emergence/decay of supersoft X-rays. 

The white dwarf in RX J0513.9— 6951 is now constantly 
growing in mass at the critical rate of Mci- We expect 
that the mass of the white dwarf increases up to near the 
Chandrasekhar mass and the white dwarf eventually ex- 
plodes as a Type la supernova. Thus, RX J0513.9— 6951 is 
the first beautiful example of the accretion wind evolution, 
which is a key evolutionary process in a recently developed 
evolutionary path to Type la supernovae (e.g., Hachisu et 
al. 1996, 1999a,b; Hachisu & Kato 2001; Langer et al. 
2000; Li & van den Heuvol 1997). 

We are indebted to F. Meyer and Emi Meyer-Hofmeister 
for stimulating us to work on RX J0513.9— 6951 during our 
stay at the Max-Planck Institute for Astrophysics. This 
research has been supported in part by a Grant-in- Aid for 
Scientific Research (11640226) from the Japan Society for 
the Promotion of Science. 
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(a) optically thick wind phase 




(b) just after the wind stops 




(c) nnoss accretion nnckes a spray 



Fig. 1. — Configurations of our modeled RX J0513.9— 6951 are illustrated: (a) in the massive wind phase (X-ray off), (b) just after the wind 
stops (X-ray on), and (c) during a rapid mass accretion phase soon after the wind stops (X-ray on). The cool component (right) is a slightly 
evolved main-sequence (MS) companion {2.6Mq) filling its inner critical Roche lobe. The north and south polar areas of the cool component 
are irradiated by the hot component (1.3 Mq white dwarf, left). The separation is a = 5.53i?0; the effective radii of the inner critical Roche 
lobes are = 1.17 Rq, and i?2 = R2 = 2.43i?0, for the primary WD and the secondary MS companion, respectively, (a) The surface of the 
accretion disk is blown in the wind and its optically thick outer edge extends up to the companion star or over. The large velocity difference 
between the wind and the disk surface certainly drives the Kclvin-Helmholtz instability at the interface, so tliat the very surface layer of the 
disk is dragged away as an almost free stream going outward. This surface free stream is initially optically thick but becomes optically thin 
outside because of geometrical dilution effect. We regard tiie transition from optically thick to thin as the edge of the disk, (b) The disk 
shrinks to a normal size (0.7 — 0.8 times the Roche lobe radius) in several dynamical timescales, i.e., several orbital periods after the wind 
stops, (c) A rapid mass accretion maJses a spray around the disk edge. 
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Flo. 2. — A modeled light curve for a limit cycle is plotted. A rapid mass accretion to the WD begins at A. The envelope of the WD 
expands to blow an optically thick wind at B. The strong wind causes a quick rise of the optical light curve, because the disk expands to a 
few to several times the previous size (see disk configuration (a) in Fig. 1). The massive wind certainly absorb supersoft X-rays. Thus, this is 
the optical high/X-ray off state. Because the wind hits and then strips-off the very surface layer of the MS companion, mass transfer from the 
MS is heavily suppressed and stops at C. However, the mass accretion to the WD still continues during a viscous timcscale of the accretion 
disk, i.e., t^is- After a large part of the mass in the accretion disk is drained, the mass accretion to the WD virtually stops at C. The mass 
of the WD envelope is now gradually decreasing due mainly to wind mass loss. Then the WD photosphere is also gradually shrinking. This 
makes a slow decline in the optical light curve. The wind eventually stops at D. The disk shrinks to a normal size at E, causing a sharp drop 
of the optical light curve (see disk configuration (b) in Fig. 1). The mass transfer from the MS resumes at E, and then the edge of the disk 
flares up at F (see disk configuration (c) in Fig. 1). It takes several dynamical timescales to change from E to F. We expect copious supersoft 
X-rays after the wind stops. This is the optical low/X-ray on state. A rapid mass accretion to the WD resumes at G(=A), because it takes 
a viscous timescale, i.e., tvis, that the matter reaches the WD surface via an accretion disk. This cycle is repeated. 
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(a) mass accretion on (b) mass accretion off 

Fig. 3. — Schematic illustrations for the flow from the secondary main-sequence companion: (a) mass transfer is going on when the stripping 
effect by the wind is smaller than the original mass transfer rate, i.e., Mgtrip < ^MS- (b) Mass transfer stops if the stripping effect by the 
wind is larger than the original mass transfer rate, i.e., Mstrip > -^MS- We are able to know when mass transfer resumes again by checking 
whether or not Mflow > by integrating equation (16), where Mflo^ indicates the level of flow (Mflow = at full). 
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Fk;. 4.-- Numerical results of our RX J0513.9— 6951 models arc shown against time (JD 2,448,000+) togctlicr with observational light 
curve: (a) Observational V-magnitudes (taken from Alcock et al. 1996). (b) Calculated 1/- magnitudes {thick solid line) together with B — V 
color {thin solid line), (c) Mass accretion rate to the WD (Mace, thick solid) and wind mass loss rate from the WD (Mwindi f^*'* solid), 
both in units of 10~ MQyr~^ . (d) Photospheric radius of the WD envelope in units of Rq {thick solid) and surface temperature of the WD 
envelope in units of eV {thin solid). The model parameters are summarized in Table 1, i.e., Mwd = 1-3 ^MS = 5.0 x 10~®Mq yr~^, 
ci = 10.0, and tvis = 20.5 days. Direct fitting to the brightness of the observational light curves indicates an apparent distance modulus of 
(m - M)v = 18.7. 
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Fig. 5. — Calculated V-magnitude light curves are plotted against time (JD 2,448,000+) for various mass transfer rates. The model 
parameters for each case are also summarized in Table 1: (a) Mms = 0.9 X IO^^Mq yr^^. It takes about 210 days for this relatively low 
mass accretion rate that the WD envelope mass reaches 6.5 x 10~'^ Mq, where the WD begins to blow a wind, from the initial envelope mass 
of 5.0X 10-''M©. The first flat part corresponds to a part of this 210 days, (b) Mms = 1-0 x lO^^'M© yr^^, (c) Mms = 2.0 x lO-^M© yr-\ 
(d) Mms = 3.0 X 10-«Mq yr-i, (e) Mms = 7.0 x IO-^Mq yr-\ (f) Mms = 10. x IO-^Mq yr"!, (g) Mms = 15. x lO-^Mo yr"!, (h) 
Mms = 20. X 1O~^M0 yr~^. Each model has M^vd = 1-3 Mq, c\ = 10.0, and tvis = 20.5 days in common. We adopt the apparent distance 
modulus of (m - M)v = 18.7. 
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Fig. 6. — Calculated y-magnitude and B — V color light curves are plotted against time (JD 2,448,000+) for our Mwd = 1-35 Mq 
models. The model parameters are summarized in Table 1: (a) Mms = 1-0 x 10~^Mq yr~i, (b) Mms = 2.0 x 10~®Mq yr~i, (c) 
Ai^MS = 5.0 X 1Q~^Mq yr~^, (d) Mms = 10.0 x 1Q~^Mq yr:~^, and (e) observational ^-magnitudes are also added (taken from Alcock 
et al. 1996). The other parameters of ci = 10.0 and tvis = 15-0 days are common among the model light curves. We adopt the apparent 
distance modulus of (m — M)y = 18.7. 
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Fig. 7. — Calculated y-magnitudc and B — V color light curves are plotted against time (JD 2,448,000+) for our My/u = 1-2 Mq models. 
The model parameters are summarized in Table 1: (a) Mms = 1-0 ^ W^^Mq yr"-*^, (b) Mms = 2.0 X 10~^Mq yr~^, (c) observational 
V-magnitudes (taken from Alcock et al. 1996), (d) A%s = 3.0 X IO^^'A/q yr^^, (e) A%s = 5.0 X IO^^A/q yr^^. The other parameters of 
ci = 10.0 and tvis = 32.0 days are common among the model light curves. We adopt the apparent distance modulus of (m — M)v = 18.7. 
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Fig. 8. — Calculated ^-magnitude and B — V color light curves arc plotted against time (JD 2,448,000+) for our My/u = 1-1 A<f0 models. 
The model parameters are summarized in Table 1: (a) Mms = 0.6 X IO^^A/q yr"-*^, (b) Mms = 1-0 x 10~^Mq yr~^, (c) observational 
V-magnitudes (taken from Alcock et al. 1996), (d) A%s = 2.0 X IO^^'A/q yr^^, (e) A%s = 3.0 X IO^^A/q yr^^. The other parameters of 
ci = 10.0 and tvis = 45.0 days are common among the model light curves. We adopt the apparent distance modulus of (m — M)v = 18.7. 
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Fig. 9.- Calculated ^-magnitude and B — V color light curves are plotted against time (JD 2,448,000+) for a different value of ci = 5.0. The 
model parameters are summarized in Table 1; (a) Mms = O.7xlO"'^M0 yr^^, (b) Mms = l.OxlO"''Af0 yr^^, (c) Mms = 2.OxlO~''M0 yr'^, 
(d) observational ^/-magnitudes (taken from Alcock et al. 1996), (c) Mms = 3-0 x 10^*'Mq yr^'-. The other parameters of M\yd = 1-3 Mq 
and tvis = 36.0 days are common among the model light curves. We adopt the apparent distance modulus of (m — M)v = 18.7. 
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Fig. 10. — Calculated y-magnitude and B — V color light curves are plotted against time (JD 2,448,000-1-) for a different value of ci = 1.5. 
The model parameters axe summarized in Table 1: (a) Mms = 1-0 x W~^Mq yr~^, (b) observational ^-magnitudes (taJjen from Alcock et 
al. 1996), (c) Mms = 2.0 x IQ-^Mq yr"!, (d) Mms = 3.0 x lO-^M© yr"!, (e) Mms = 5.0 x lO-^M© yr"!. The other parameters of 
A^WD = 1.3 Mq and tvis = 62.0 days are common among the model light curves. We adopt the apparent distance modulus of (rn—M)v = 18.7. 
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Fk;. 11.-- Three times the Kelvin- Helmholtz timescale of the WD envelope is plotted against the WD mass. Since the timcscale depends 
on the metallieity (Z) and the hydrogen eontent (X), we plot four cases of the hydrogen content, X = 0.1. 0..3."), O.-S. and 0.7, number of 
which is attached to each line, each for two cases of the metallieity, Z = 0.004 (solid line) and Z = 0.02 (daslied- dotted line). Here we have 
calculated tkh =(thermal energy of the WD envelope)/ (luminosity of the WD) just after the wind stops. This corresponds to the shortest 
timescale that the WD photosphere shrinks by a factor of four, i.e., from ~ 0.08 Rq (~ 25 eV) to ~ 0.02 Rq (~ 50 eV). 
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Fig. 12. — Calculated y-magnitude and B — V color light curves are plotted against time (JD 2,448,000+) for various inclination angles of 
(a) i = 0°, (b) i = 20°, (c) MACHO observation taken from Alcock et al. (1996), (d) i = 28°, and (d) i = 35°. The other model parameters 
are the same as those in Fig. 4. 
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Fig. 13. — Calculated ^-magnitudes are plotted against two orbital phases (from —1.0 to 1.0) for (a) the optical high state and (b) the 
optical low state. We adopt the inclination angle of i = 28° . The other model parameters are the same as those in Fig. 4. 



